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ABSTRACT 

We give a mean-field calculation for the odd-resonating-valence-bond ORVB pairing 
scheme. We obtain interesting quasi-particle excitation energy as a function of momentum 
k. It is distinctively different from those of the d x 2_ y 2-wave, the anisotropic-s-wave, and the 
p-wave. It is a gapless theory for superconductivity with well defined Fermi surface. The 
ground state of the ORVB scheme is not an eigenstate of the parity or the time-reversal 
transformation, thus both symmetries are violated. Some of them are already manifested in 

^ -E-k- It is interesting to find out if such pairing order-parameter scheme exits in some 
materials in nature. 



The discovery of high T c superconductivity^ led to an enormous interest in models 

[21 [3—51 

with strong correlations, especially the Hubbard model with strong repulsion. In a 
previous paper , we pointed out the existence of extended local symmetries in the large-C/ 
Hubbard model at half-filling as well as the anti-ferromagnetic exchange model, and proposed 
a new order-parameter scheme, the odd-resonating-valence-bond (ORVB) scheme. We later 
in Ref.[7] showed that, besides the current popular d x 2_ y 2-wave order-parameter scheme, the 
ORVB scheme is another natural consistent solution to the rigorous constraints imposed on 
the system under the thermal equilibrium. (Earlier Zhang pointed out that the pure s-wave 
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order-parameter scheme is not allowed. ) 

In this paper, we report the result of a mean-field calculation for the ORVB pairing 
scheme. We obtain interesting quasi-particle excitation energy as a function of the 
momentum k. It is a gapless theory for superconductivity with well defined Fermi surface. 
They are distinctively different from those of the d x 2_ y 2-wme pairing, the anisotropic-s-wave 
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pairing, the p-wave pairing, and many other pairing schemes. Our ORVB scheme can 
be tested in the angle-resolved-photo-emission-spectroscopy (ARPES) experiments, however 
the experiments must scan many key directions in the whole Brillouin zone. The ground 
state of the ORVB scheme is not an eigenstate of the parity or the time-reversal operator, 
thus both symmetries are violated. Some of these symmetry violation effects have already 
manifested in ^ E_^. 

It is interesting to find out experimentally if such pairing order-parameter scheme exits 
in some materials in nature or if it is a scheme for some high T c superconductivity . 

We start with the following mean-field Hamiltonian of the tight-binding model with 
anti-ferromagnetic exchange, 

Hmean = ^(r,r'),a(cJ j0 -C r ' jf7 + h.C.) + //E r)CT cJ ;f7 C r;CT 

" (<//2)£( r ,r0{Ar,r'Aj. r , + A* r ,A r , r , + X r,r>xl,r> + Xr,r'Xr,r'}, (1) 



where 

A r ,r' = Cr.TCr'.J. ~ Cr,|Cr',T , Xr,r> = 4,T £ r',T + 4,| £r U ■ ( 2 ) 

These two operators are related by the SU(2) symmetry, which is the symmetry of the 
J-term or the symmetry of the anti-ferromagnetic exchange. In obtaining the mean-field 
Hamiltonian, we have used the ORVB order parameters: A r>r / =< A r>r / >= Tr(e^ H A rjV >), 
where H is the Hamiltonian and (3 = (kT)^ 1 . At temperature T = 0, < A r>r ' >=< 
0|A r]I .'|0 >, i.e., the expectation value of the ground state; similarly for Xr,r' =< Xr,r' >■ 
The doping parameter 5 is defined by 1 — 5 = N~ 2 S r>cr < c\^c r ^ a > . 

The ORVB order-parameters can be parameterized in the following way: 

A r , r > =V,r±e e™ Ae ie , (3) 

Xr,r' =V,r±e X***, (4) 

where the unit vector e = e x or e y . As shown in Ref. [7], the rigorous constrains imposed 
by thermal equilibrium allow all these real order parameters, A, 9, x and <fi to be nonzero 
and unequal in different directions. For simplicity, as done in many other order parameter 
schemes, we choose them to be constants, independent of sites r. Notice that from the 
definitions of A r>r ' and Xr,r' in Eq.(2), we have A_ r _ r ' = — A r r ' and x-r-r 1 = (Xr,r')* • In 
the case A = but x ^ and <fi ^ 0, the ORVB scheme becomes the flux phase scheme! 10 ^ 
However the main new features of the ORVB scheme come from A ^ and the interplay 
between the two sets of order parameters given in Eqs. (3) and (4). 

After the Fourier transformation, the Hamiltonian in momentum space becomes 

H°™ n B =Z Ka {2tC k +11- J X [cOs(k x + 0) + COs(k y + 0)]}4,a £ k,a 

- £ k {JA^e-^ + %_ k , T c ka + e^ci/,^]}, (5) 
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where C k = cosk x + cosk y and S k = sink x + sink y . Obviously from Eq. (5) we see that the 
momenta of the electrons in the pairing are k and tx — k, in contrast to the k and — k of the 
BCS pairing! 111 

The Hamiltonian Hmean B can be diagonalized by the following Bogoliubov transforma- 
tion c k>i = u* k 7k,| + v k 7t_ k T and ct_ k;T = ~v* k 7k, j + u k , where % a and 7^ are 
quasi-particle operators. Then the Hamiltonian becomes 

frORVB _ v (rpORVB-i ~ , pi(h / r \ 

with the energy given by 

E ORVB =(2t _ j XCOS(f)) c k ± ^/[JAS fc ] 2 + [// + J*Sfcsm0] 2 , (7) 

where ± are taken as the sign of [// + JxSksin(f>}. Notice that when the order parameters 
become zero, E^ RVB = e k = 2tC k + [i ■ 

The ORVB energy expression, Eq. (7), is quite different from those of the BCS-type 
of theories, E BCS = ±A/e k 2 + [JA k ] 2 , in which A k = A s for the s-wave scheme (the 
original BCS order parameter); A k = A^(cosk x — cosk y ) for the d x 2_ y 2-wsve scheme; 
A k = A a [(cosk x — cosk y ) 4: + constant] for the anisotropic s-wave scheme; (the sign in front of 
fiBCS j g f-^at of e k ). Notice that when the order parameters are zero (at T > T c ), the ORVB 
scheme has the same energy expression, E k = e k , as the BCS-type of theories. When the 
order parameters becomes nonzero (at T < T c ), the BCS-type of theories gives energy discon- 
tinuity (gap) around e k = 0. The ORVB scheme gives very different energy discontinuities, 
which we shall discuss later. 

Notice that the momentum-factor associated with A in the E^ RVB of Eq.(7) is S k , 
a result of the pairing momentum being ir. The momentum-factor associated with x is 
cos(k x + 0) + cos(k y + 0). It decomposes into CkCos<f> and Sksincj); the former gives the term 
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outside the square root in Eq.(7), and the latter gives the term inside the square root in 
Eq.(7). 

The ground state is defined by 7 k;(T |0 >= for > and 7 kf7 |0 >= for < 0. 
The quasi-particle operators 7 k g and 7 kiCr are free fermions. < 7 k a 7k,o- > are the only non- 
vanishing expectation- values. The requirement that the occupation number at a temperature 
T is given by the Fermi-Dirac distributions < 7 k(J 7k,o- >= [e^ £k + imposes a set of 
consistent conditions upon the pairing parameters. After Fourier-transforming to momentum 
space and then Bogoliubov-transforming c and c 1 to 7 and , we obtain the following four 
equations: 

5 = N~ 2 S k [J X S k sin(f> + ^X' 1 , (8) 



where 



2J X 2 cos 2 <j ) - At X cos<t> = N~ 2 S k (F k G k ), (9) 

2J (A 2 + X 2 sm 2 0) + ^ = N' 2 S k (X k F k ), (10) 
2 J(A 2 - x 2 sm 2 0) - fi5 = 

N- 2 S k {[JA,S fc ] 2 - [JxS k sin<f> + /i] 2 }F k X k \ (11) 



X k = ^[JAStf + lJxSkStn^ + fi] 2 , (12a) 

F k = ( J X cos<t> - 2tS)C k: (126) 

F k = sinh(pXk)[cosh(pX±) + cos/i(/3y k )] _1 , (12c) 

G k = sin/i(/?y k )[cos/i(/3X k ) + cosh^Y^)]' 1 . (12d) 



For given temperature /3 and doping parameter 0", totally there are four equations, 
Eqs.(8), (9), (10), and (11), for the four unknowns A, x, 0, and (Notice that the phase 6 
does not appear in these equations.) 



At T = 0, the temperature dependent factor F k — 1 for X k > |Y k | and F k = for 
Xk < \Y k \; G k = 1 for F k > X k , G k = -1 for F k < -X k , and G k = for |y k | < X k . 
Therefore Eqs.(8), (9), (10), and (11) are further simplified, which can be trivially obtained 
and we do not write them out explicitly. 

We have analyzed Eqs.(8), (9), (10), and (11) numerically at small values of the doping 
parameter 5 at various values of T. For a given small value of 5, we find that the solution for 
these equations giving a small value of ji and Eq.(ll) does not impose a strong independent 
condition from that of Eq.(10). [These can be understood by the fact that at 5 = and 
[i = Eq.(ll) reduces to Eq.(10) and Eq.(8) becomes an identity. However, our solution is 
different from the one exactly at half-filling, i.e. 5 = and fi = 0, where further degeneracy 
happens and gives a different solution.] Therefore for a small doping parameter S only two of 
the parameters, x 2 cos 2 and A 2 + x 2 sin 2 (f) are determined, essentially by Eq.(9) and Eq.(10). 
The approximate solutions we have obtained are the following: the critical temperature T c 
(that we define to be the temperature below which all the order parameters are nonzero) 
and the order parameters at T = are all of the order of the interaction strength J; the 
order-parameters decrease as the doping parameter 5 and/or as the temperature T increases; 
T c decreases with the increasing value of the doping parameter 5. For example, at T = we 
find the approximate solution at 5 ~ 0.18 (t ~ 0.18): /i ~ 0.08, A 2 + x 2 sin 2 4> ~ 0.032, and 
X 2 cos 2 (f) ~ 0. From these discussions, we also see that our ORVB scheme contains the flux 
phase (i.e. x and A = 0) and the case of x — and A ^ as special cases. We shall 
see below that it is the general cases of both x 7^ and A^0 that give the interesting new 
physical phenomena. 

With these solutions, we can calculate the quasi-particle excitation energy E^ RVB of 
Eqs.(6) and (7). In Fig.l we show the E^ RVB distribution as function of k for A ~ 0.17 
and — —0.06. Notice that the lattice lines of the Brillouin zones are no more the 

symmetry lines. We shall show later, this is a result from the violation of the parity as well 
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as the time reversal symmetries of our ORVB pairing scheme. The xsm0 7^ gives a wall of 

discontinuity (gap) in the k-plane specified by Sk = —fj,/(Jxsin<f>), surrounding the centers 

at (k x , k y ) = (vr/2, vr/2), modulus 2n in both k x and k y . The magnitude of the discontinuity 

[121 

is \(2fiA)/(xsin(f))\. Notice in Fig. 1 there are two lines of zero energy (gaplessness). The 
line of -Ek = ek is also symmetrically centered around (k x ,k y ) = (n/2,n/2), which is not 
shown in the figure but can be easily seen from Eq.(7). 

The quasi-particle excitation energy Ek distributions for the BCS-type theories are very 
different from that of the ORVB scheme we just presented. For BCS-type theories, the 
lattice lines of the Brillouin zones are the symmetry lines, a result of parity and time reversal 
symmetries. The anisotropic s-wave scheme is a finite-gap theory (i.e. E^ 7^ in the whole 
Brillouin zone) and the d x 2_ y 2-wave scheme has isolated four points of gaplessness. We shall 
mention at the end of the paper how these differences can be experimentally distinguished. 

At T = 0, the occupation number in the quasi-particle space is = < 7^ |7k,| + 1^ 1 7k, j, > 
= for E^ > and = 2 for E^ < . After Bogoliubov-transforming 7 and 7^ back 
to c and c\ we obtain the occupation number in the electrons n| =< cj. |Ck,| + , Ck,j, >, 
shown in Fig. 2. These occupation number distributions are also very different from those of 
the BCS-type of theories. 

We note that for every %sin0 < solution, %sin0 = —0.06 in Fig.l, there is another 
equally allowed ORVB solution with xsin§ > with the discontinuity center changed to 
(k x , k y ) = (—it/2, — 7r/2). The other characteristic are the same as those just discussed and 
presented for %sin0 < 0. 

Besides the ORVB, there is another kind of solution to the thermal equilibrium constraint 
equations, which we call the (ORVB)' pairing scheme. In stead of Eqs.(3) and (4), we have 

A r , r , =V,r+e x e^A e w = -V,r+e y A e id , (13) 
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Xr,r' =V,r+e x X ^ = Sr> ,r+e y X ^ , (14) 

and A, x, 9, and <fi are real constants. The net result is that for xsinfy < 0, the discontinuity 
shifted to surround the centers at (k x ,k y ) = (n/2, —ir/2) modulus 2n in both k x and k y ; 
for xsin§ > 0, the discontinuity shifted to surround the centers at (k x ,k y ) = (— 7r/2,7r/2) 
modulus 27r in both k x and k y . 

Next, we discuss the intrinsic parity and time reversal symmetries. Both the A r)I .' and 
Xr,r', Eqs.(3) and (4), violate the parity and the time reversal symmetries. 

Under parity transformation P, for any eigenstates of parity \a > and \b >, the expecta- 
tion value of any operator O has the identity < a\0\b >=< a|POP _1 |6 > . From the fact 
Pc r , a P~ l = C- r ,a and Pcl^P^ 1 = cL r;Cr , we obtain Ph. r ^P^ 1 = A_ r; _ r / and PXr,v'P~ l = 
X-r,-r' from the definitions of these operators in Eq.(4). If the eigenstates of Hamiltonian 
are eigenstates of the parity operator, the first relation implies < A r r / >=< A_ r _ r / >; in 
which the right-hand-side (r.h.s.) is A r;I .', and the left-hand-side (l.h.s.) is — A rjI .', using 
Eqs.(3) and (4); thus must A r)1 ./ = 0. The second relation implies < x>,r' >=< X-r-r 1 >, in 
which the r.h.s. is Xr,r' and the l.h.s. is xj? r 'j thus must Xr,r' = real, i.e., = 0. Therefore, 
A 7^ or the complexity of x e ^ \ i- e - ; 4> 7^ ; implies that eigenstates of the Hamiltonian 
can not be the eigenstates of parity, i.e., parity is violated. Notice that this parity violation 

property is different from that of the p-wave pairing scheme, in which the ground state is an 

[9] 

odd eigenstate of the parity transformation. 

Under time reversal transformation T, T\a >= \a t > and T\b >= \b t >, the expectation 
value of any operator O has the identity < a\0\b >=< at\TOT~ l \b t >* . From the fact that 
Tcr^T -1 = c r ^ a and Tcl^T -1 = cl_ a , we obtain TA r)r /T _1 = — A r r / and Txr,v'T^ 1 = 
Xr,r' from the definitions of these operators in Eq.(4). If the ground state is the eigenstate of 
time reversal T|0 >= |0 >, < 0|A r;I ./|0 >=< 0| — A r?r / jO >* and < 0|xr,r']0 >=< 0|xr,r'|° >*; 
which implies at temperature T = 0, A rir ' = — A* r , (i.e., A rjl .' is pure imaginary and 
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9 = 7r/2) and Xr,r' = Xrr' (i- e -> Xr,r' is real and = 0, n). Thus AcosO ^ or xsin(j) ^ 
implies the violation of time reversal symmetry. (Note that — Xi being real, is not 
an interesting order parameter, Since it only renormalizes the tight-binding term in the 
Hamiltonian.) 

If the ground state is an eigenstate of time reversal T|0 >= |0 >, following similar 
discussions as above, we can easily show that =< 0|ff mean |0 >= £Lk from either the 
Hermiticity of the Hamiltonian, H mean = Hmean, or TE mean f^ = H mean . Therefore, £ k ^ 
indicates that the ground state is not an eigenstate of time reversal. In our case, 
£"k ^ comes solely from the xsin^ term. If %sin0 = 0, our E^ = E_^ even if there is 
time-reversal violation from AcosO ^ 0; thus the time reversal violation effect from AcosO 
can not be observed in E^, since E^ is not sensitive to the phase of Ae t0 . 

Finally we point out some distinct features of the ORVB scheme that experiments can 
look for. In our ORVB scheme the wall of the energy discontinuity is off-centered in the 
Brillouin zone, centered in one of the quadrant of the Brillouin zone (see Fig. 1 and discus- 
sions in the text around the Eqs.(12) to (14)), in contrast to those of the cL^^-wave and 
the anisotropic s-wave which are all centered symmetrically around the center of the Bril- 
louin zone. It will be interesting to study such asymmetric distributions in experiments that 

[131 

can measure the energy discontinuities, like the ARPES experiments. Due to complexity 
of the order parameters, in experiments that can measure the phases of order parameters 
(like the Josephson interference experiments) the ORVB scheme can give interference phase 
shifts different from 7r, which is the only kind of interference phase shift the <i a .2„ y 2-wave 
scheme can produce (due to the reality of its order parameter). We have also calculated the 
temperature dependence of the penetration depth. It increases with temperature faster than 
that from the d^^-wave case. (We shall publish our detail numerical results in a separate 
paper.) 
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To conclude, our ORVB scheme offers a pairing order parameter scheme for supercon- 
ductivity very different from the current popular models, e.g., the d x 2_ y 2-wave scheme, the 
anisotropic s-wave scheme, and the p-wave scheme. It has interesting and distinct features 
that can be experimentally tested. Seeing the ORVB characteristics in some materials will 
certainly be very exciting. 
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FIGURE CAPTIONS 

Fig.l The equal E^, the quasi-particle energy, contour diagram as function of the momen- 
tum k in the ORVB scheme with a small doping parameter 5. For the values of parameters 
used and discussions, see the text. Notice that the wall of discontinuity (gap) is centered 
around the center of a quadrant of the Brillouin zone. Shown in the figure is the case of 
the center being at (k x ,k y ) = (tt/2, tt/2). It is equally likely in the ORVB scheme that the 
center is at (-tt/2, tt/2), (tt/2, -tt/2), or (-tt/2, -tt/2). 

Fig. 2 The three dimensional view of the electron density n| of the ORVB scheme. 



10 



REFERENCES 

[I] J. G. Bednorz and K. A. Muller, Z. phys. B64 (1986) 189; M. K. Wu, J. R. Ashburn, 
Y. Q. Wang, and C. W. Chu, Phys. Rev. Lett. 58 (1987) 908. 

[2] J. Hubbard, Proc. Roy. Soc. A276 (1963) 238; P. W. Anderson, Science 235 (1987) 
1196. 

[3] G. Baskaran, Z. Zou, and P. W. Anderson, Solid State Commun. 63 (1987) 973; V. 
Emery, Phys. Rev. Lett. 58 (1987) 2794; G. Kotliar, Phys. Rev. B37 (1988) 3664. 

[4] S. Chakravarty, A. Sudb0, P. W. Anderson, and S. Strong, Science 261 (1993) 337. 

[5] P. Monthoux and D. Pines, Phys. Rev. B47 (1993) 6069; D. Pines, Physica B199 and 200 
(1994) 300; D. J. Scalapino, Physics Reports 250(1995) 329; J. R. Schrieffer, Sol. State 
Comm. 92 (1994) 129. 

[6] L.-L. Chau, D.-W. Huang, and Y. Yu, Phys. Rev. Lett. 68 (1992) 2539. 

[7] L.-L. Chau and D.-W. Huang, Constraints on Order Parameters and Correlation Func- 
tions of Systems in Thermal Equilibrium, submitted to Phys. Rev. B. 

[8] S. C. Zhang, Phys. Rev. B42 (1990) 1012. 

[9] For example the OPSP pairing scheme of R. T. Scalettar, R. R. P. Singh, and S. C. 
Zhang, Phys. Rev. Lett. 67 (1991) 370; the p-wave pairing scheme of A. V. Balatsky, E. 
Abraham, D. J. Scalapino, and J. R. Schrieffer, Physica B (1994) 363. 

[10] I. Affleck and J. B. Marston, Phys. Rev. B37 (1988) 3774. 

[II] J. Bardeen, L.N. Cooper, and J.R. Schrieffer, Phys. Rev. 108 (1957) 1175. 

[12] As xsin$ — > 0, the discontinuity disappears and the energy E^ becomes continuous: 
j^ORVB _ — Jxcos(f))Ck + ^/[JASk] 2 + At 2 . As \xsin<f)\ increases (remember that A 2 + 
X 2 sin 2 (fi ~ constant), the area enclosed by the discontinuity enlarges and the discontinuity 

11 



along the boundary becomes smaller. When \xsiri(p\ reaches its maximum (A = 0), the 
spectrum becomes continuous again, E^ RVB = (2t — Jxcos<f))Ck + [i + JxSk^ 11 ^- 

[13] It is impressive that the energy discontinuities can be measured at all. So far only 
partial regions of the Brillouin zone of Bi2212 has been scanned for the energy discontinuities 
by the ARPES experiments, see Z.-X. Shen, et al, Science 267 (1995) 343, and references 
therein. 



qc;manuscripts;orvb; O_orvb_tb_cond.tex 11c :6/30/95 



12 



